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Abstract 


We analyzed the spectral properties and pulse profile of PSR J1811-—1925, a pulsar located in the center of 
composite supernova remnant (SNR) G11.2-0.3, by using high timing resolution archival data from the Nuclear 
Spectroscopic Telescope Array Mission (NuSTAR). Analysis of archival Chandra data over different regions rules 
out the SNR shell as the site of the hard X-ray emission while spectral analysis indicates that the NuSTAR photons 
originate in the pulsar and its nebula. The pulse profile exhibits a broad single peak up to 35 keV. The jointed 
spectrum by combining NuSTAR and Chandra can be well fitted by a power-law model with a photon index of 
T= 1.58 +0.04. The integrated flux of jointed spectrum over 1-10keV is 3.36 x 107" erg cms a Aa 
spectrum of pulsar having photon index T = 1.33 + 0.06 and a 1-10 keV flux of 0.91 x 107'? erg cm7 ' We 
also performed the phase-resolved spectral analysis by splitting the whole pulse-on phase into five phase bins The 
photon indices of the bins are all around 1.4, indicating that the photon index does not evolve with the phase. 
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1. Introduction 


PSR J1811-1925 is an X-ray pulsar located at the center of 
supernova remnant G11.2—0.3, a young core-collapse super- 
nova remnant (CC SNR) (Madsen et al. 2020). The system is 
composed of a shell emitting soft X-ray radiation, a pulsar at 
the center with pulsar wind nebulae (PWN) around (Dean et al. 
2008). In the radio image by VLA and X-ray image by 
Chandra, the outer boundary of G11.2-0.3 shows a near- 
circular structure (Borkowski et al. 2016). G11.2—0.3 is one of 
the youngest CCSNRs in the Milky Way with age ~3600 yr 
(Torii et al. 1999). The supernova remnant is associated with 
the “guest star” event recorded in AD 386 (Clark 1978). 
However, the high visual extinction toward G11.2-0.3, 
obtained both from [Fe II] observations and from absorption 
measurements toward the PWN and the shell, rules out the 
association of G11.2-0.3 with any naked-eye event seen on 
Earth (Koo & Lee 2015; Borkowski et al. 2016), in particular 
with the AD 386 event. According to the measurement of 
nebular expansion rate in the X-ray range by Chandra 
Observatory, the age of the SNR is estimated to be between 
960 and 3600 yr (Tam & Roberts 2003). The out-flowing wind 
is non-spherical and more likely to be limited to a bipolar 
outflow (Roberts et al. 2003). 

PSR J1811—1925 was discovered in X-rays by ASCA, with 
spin period P = 65 ms and Ê = 4.4 x 107!4 s s~!. Its magnetic 
field is about B ~ 10'? G, and its spin-down energy loss rates 
È ~ 6 x 10° erg s~! (Torii et al. 1997, 1999). The pulse signal 
is detected both in near-infrared (Moon et al. 2009) and X-ray 
band (Kaspi et al. 2001). But no pulse signal was detected in GeV 


by Fermi-LAT (Acero et al. 2016) and in the radio band (Tam & 
Roberts 2003). The characteristic age of PSR J1811-1925 is 
23.3 kyr, which is significantly older than G11.2-0.3 (Kaspi et al. 
2001). The flux is (1.22+0.21) x 107''ergcem™’s~' in the 
soft y-ray (20-100 keV) band by INTEGRAL/ISGRI (Kuiper 
& Hermsen 2015). It is suggested that the hard X-ray radiation 
comes from the pulsar in the center of the supernova remnant. 
Using the high sensitivity of the Nuclear Spectroscopic 
Telescope Array Mission (NuSTAR) in the hard X-ray energy 
band (3-79 keV), we can explore the radiation mechanisms of 
pulsars in different emitting regions in more detail. In addition, 
the phase-resolved spectroscopy of the pulsar is a crucial 
method to study the emission generation mechanisms, and the 
phase-resolve spectrum of PSR J1811-1925 has not been 
reviewed yet. PSR J1811-1925 is a good sample for 
understanding the soft y-ray pulsar population. 

In this work, we introduced details of data reduction in 
Section 2. By combining data from NuSTAR and Chandra, we 
obtained the pulse profile, the jointed spectrum of PSR 
J1811-1925 and its PWN, pulse spectrum and phase-resolved 
spectroscopy of PSR J1811-—1925, and these are presented in 
Section 3. Finally, we discuss and interpret these results in 
Section 4. 


2. Observations and Data Reduction 
2.1. NuSTAR Observation and Data Reduction 


NuSTAR is a hard X-ray telescope with two instruments 
(usually labeled by their focal plane modules, FPMA and 
FPMB), operates in the band from 3 to 79 keV with excellent 
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Figure 1. NuSTAR FPMA 3-79 keV image of the region around PSR J1811—1925. Black circles: (a) radius of circle is 60”, where the spectrum of source extracted; 
(b) radius of circle is 74”, the spectrum of background extracted from the region between (a) and (b). 


Table 1 
Information of Observations 
Name ObsID Date MJD Exposure Time Instrument 
G11d2m0d3 90201027002 2016-06-22 22:11:08 57561.92 90 ks NuSTAR 
SNR GO11.2-00.3 14831 2013-05-05 19:28:01 56417.81 180 ks Chandra/ACIS-S 


spectral resolutions (400eV @ 6keV, 900eV @ 60 keV) 
(Harrison et al. 2013). PSR J1811-1925 was observed by 
NuSTAR around June 22 2016 (OBSID 90201027002) with an 
exposure time of 90 ks, see Table | for more details. 

In this work, we analyzed NuSTAR data using HEASoft 
(version 6.28), standard nupipeline and nuproducts 
utilities (Version: 0.4.7) and NuSTAR CALDB v20210104. 
We extract the source spectrum from an circular region of 
radius 60” centered on the pulsar position at a = 18"11™29822 
and 6 =— 19252756 by DS9 (version 4.1) (as shown in 
Figure 1), the background spectrum is extracted from a annular 
region of 60”—74” (Kaspi et al. 2001). The arrival time of each 
event is corrected to solar system Barycenter (SSB) by tools 
barycorr with solar system ephemeris DE200. We extracted 
phase-resolved spectroscopy by dividing the whole pulse-on 
phase into five phase bins by using XSELECT, the background 
spectrum is extracted at phase regions (phase 0.6—0.7). The 


pulse spectrum of the pulsar is extracted at phase bins with 
prominent pulse signal. 


2.2. Chandra Observations and Data Reduction 


Chandra observed SNR G11.2-0.3 12 times from 2000 to 
2013, with a wide field of view covering the remnant SNR 
G11.2-0.3. (Borkowski et al. 2016). We selected data of ObsID 
14831 based on its long exposure of about 180 ks and 
observation time close to that of our NuSTAR data (see Table 1 
for more details). In the observation of ObsID 14831, the pulsar 
is positioned on the chip S3 of the Advanced CCD Imaging 
Spectrometer (ACIS) in the standard exposure mode (Kaspi 
et al. 2001). 

The Chandra data are analyzed and calibrated using the 
Chandra Interactive Analysis of Observations (CIAO V4.14) 
software package and the Chandra Calibration Database 
(CALDB V4.9.8), respectively. PSR J1811-1925 is a bright 
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Figure 2. Chandra image of the region around PSR J1811—1925. White circles: (a) 03-60”, where the source spectrum is extracted; (b) 60”—74", where the 
background spectrum is extracted. The spectra of pulsar and its background are extracted from annular region (region (c) : 03—10) and annular region (region (d) : 


10-138) at the center region. 


point source. To avoid the pile-up effect, we extracted the 
jointed spectrum of the source and background from two 
separate annular regions: region (a) with a range of 0” 3-60” (as 
shown in Figure 2), and region (b) with a range of 60”—74”. See 
Figure 2 for visual reference. With an excellent angular 
resolution (about 0”5) of Chandra, the pulsar can be resolved 
from Chandra image. The spectrum of the pulsar can be 
extracted from an annular region (0”3-—1”0, region(c)), and the 
spectrum of background is extracted from an annular region 
(170-1138, region (d)) (as shown in Figure 2), the selections of 
the regions are same with Roberts et al. (2003). 


2.3. Timing Analysis 


The time resolution of our data from Chandra observation is 
about 3.2s, which is not suitable for timing analysis (Kaspi 
et al. 2001). Here, we only perform timing analysis using 
NuSTAR data. 

In our timing analysis, to maximize the signal-to-noise ratio 
of pulse profile, all events of FPMA and FPMB in 3-79 keV 
are selected to generate a integrated pulse profile. All events 
within source region “(a)” in Figure 1 are used for spin 
frequency searching and time of arrivals (TOA) of all these 
events are corrected to SSB. The pulse phase @ of each events 


arrived at time ¢ is calculated by using a Ist-order Taylor 
expansion: 


$ = po + V(t — to), (1) 


where ġo and v are phase and spin frequency at the reference 
time fọ. Then we obtained a pulse profile using all these 
photons with their phases. The best spin frequency is a 
frequency that makes the folded pulse profile deviate the most 
from a uniform distribution as represented by the Pearson x? 
test. We search the best spin frequency by using this method 
and obtained a best spin frequency of v = 15.4564269(1) Hz at 
the epoch of To = 57563.0000001 MJD (Modified Julian Date), 
which is consistent with the period measured by Madsen et al. 
(2020). We plot our pulse profile in Figure 3. 


2.4. Spectral Fitting 


The phase-averaged spectra are generated from Chandra and 
NuSTAR data. For NuSTAR data, we use grppha from FTOOL 
in order to obtain a better signal to noise ratio (S/N) spectrum. 
Both spectra of FRMA and FPMB are rebinned using grppha to 
a minimal bin of 16, 32, 128, 256 for channels range 1-224, 
225-480, 481-1248 and 1249-2016, respectively. For these 
spectra, we use x? Statistics for spectral fitting. For Chandra 
data, source and background are extracted from nearly the 
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Figure 3. X-ray pulse profiles of PSR J1811—1925 in six regions, two pulse periods are presented for clarity. 


identical region as NuSTAR data (region (a) and (b) in Figure 
2). The spectrum of the source from Chandra is rebinned using 
grppha to a minimal bin of 3 and 20 for channels range 
69-137 and 138-697, respectively. The same group parameters 
are also applied to the pulse spectra of the pulsar and the phase- 
resolved spectra. The phase averaged spectrum is fitted with a 
Tuebingen-Boulder ISM absorption model multiplied a power- 
law plus a plane parallel shock model (see Borkowski et al. 
2001 for details of this model), given in XSPEC as: tbabs 
(powerlaw+vpshock) (Borkowski et al. 2001). We produce the 
spectrum of the pulsar by combining Chandra (from an annular 
region of 073-1”0) and NuSTAR data (from the pulse-on 
phase). The phase-resolved spectra and the spectra of pulsars 
are fitted with a simple tbabs (powerlaw) model. 

All the spectra are analyzed with XSPEC (version 12.11.1) 
(Arnaud 1996). Unless otherwise noted, the errors of spectral 
parameters are at 90% confidence level. 


3. Results 
3.1. Pulse Profile 


With our best spin frequency of v = 15.4564269(1) Hz, we 
generated six pulse profiles by extracting photons from six 
different regions, 0”—37”, 37’-60", 60”—90”", 90”—123”, 123”— 
150”, 150”-180” from NuSTAR observations, to check 
whether the pulse signal exists. As shown in Figure 3, there 
is no apparent pulsed signal beyond a radius of 150”, which is 
consistent with result of Roberts et al. (2003). Figure 4 presents 
the pulse profiles over four energy ranges 3.0-6.0, 6.0—-10.0, 
10.0-35.0, and 35.0-79.0keV. The phases of single pulse 
peaks remain approximately constant throughout the four 
energy ranges, and all normalized pulse profiles show a single 
peak that is stable except for an increase in noise observed at 
the pulse profile at 35-79 keV. 


3.2. Spectral Analysis 
3.2.1. Whole Spectrum Around PSR J1811—1925 


The whole spectra obtained from Chandra and NuSTAR 
could be roughly fitted by tbabs (powerlaw+vpshock). As 
listed in Table 2, the ionization timescale of vpshock model 
is frozen to be T = 4.2 x 10'' scm ° (Roberts et al. 2003). The 
Miota index of the spectrum is 1.58 + 0.04 with the reduced 

Ka Xa dof = 1.21 1/ 80. From the D model, the absorbed 


flux is 3.3670 x 107!2 erg cm~? s7! in 1-10 keV, which is 
consistent with the result of Gavriil et al. (2004). The emission 
lines of Mg, Si could be seen clearly from the spectrum, while 
no evidence for S line in this spectrum. The missing of S line 
might be due to the fact that the area we selected is different 
from previous studies. We found the existence of S line by 
generating spectra in different regions, and this result is 
consistent with Roberts et al. (2003), in order to ensure that 
they were as consistent as possible with the region of NuSTAR, 
the S line was not detected clearly. The presence of element 
lines means that the spectrum contained the contribution of 
SNR. Panel (b) of Figure 5 shows the ratio between the model 
and the data. From the ratio, the Mg, Si lines are not well-fitted 
enough. Therefore, we do not analyze the individual abundance 
values quantitatively and make general comments only instead. 


3.2.2. Spectrum of Pulsar 


Considering the high angular resolution of Chandra, the 
spectrum of pulsar is generated using region (c) in Figure 2, 
which is dominated by the central pulsar, the background 
spectrum is generated using region (d) around the region (c) (see 
Figure 2). For NuSTAR data, to minimize the contamination of 
SNR, we generated the spectrum of pulsar by using photos at the 
pulse-on phase, and generated background spectrum using 
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Figure 4. X-ray pulse profiles of PSR J1811-1925 in four energy ranges. We show two pulse periods for clarity, and pulse emission has been detected up to 35 keV. 


3.2 x 10°’ atoms cm’. The photon index is 1.33 +0.06 with 
the reduced x 1.14 (d.o.f 70). Cremar is the ratio different 
between the spectra of NuSTAR and Chandra and its value is 


Table 2 
The Spectral Fitting Results of Phase Averaged Spectrum by Combining 
Spectra of Chandra and NuSTAR 


No. Parameter Value 1.67, which is slightly large as the method of generating 

1 Ng? 2.50 + 0.28 spectrum of Chandra is different from NuSTAR’s as described 

2 T 1.58 + 0.04 in methods. From the power-law model, the absorbed flux is 

4 Norm (10$ 7.222037 (0.91 +0.04) x 10`! ergcm™?s™! in 1-10keV, which is 

5 kT ne 0.21704 consistent with the result of Roberts et al. (2003). 

6 Mg 0.49939 

7 si? 1.11737 

: EA me 3.2.3. Phase-resolved Spectroscopy 

6 Flux“ (1-10 keV) 3.367595 With the high timing resolution of NuSTAR, we also studied 

T Cremas 0.94 +0.04 the phase-resolved spectroscopy by dividing the whole period 

8 Xrea(d-0-f) 1.21 (80) into ten phase bins. The spectrum of the pulse-off phase 
(0.4-0.9 in Figure 8) are combined to be one spectrum because 

Notes: the counts rates of a single phase of them is too low to obtain 


è Unit: 10°? atoms cm77. 
P Unit: Solar Abundance. 
3 


enough statistics. So we obtained six spectra finally for 
studying phase-resolved 


spectroscopy. Spectra from the 


. . 11 =3 
€ Parameter is fixed, unit: 10'' scm”. 
2-1 


d The absorbed X-ray flux unit is 107}? erg cm™? s different phases were grouped with a minimum of 50 counts 


per spectral channel and fitted using XSPEC. We fitted all six 
spectra using a tbabs (power law) model with Ny fixed at 
3.2 x 10”? atoms cm~*, since NuSTAR has poor restrictions on 
the Ny value. The fitted results are listed in Table 4, and the 
residuals of the spectra fitting are plotted in Figure 7. From 
Table 4, X, 4 of all spectra show that the model can fit the 


photons at the pulse-off phase. The pulsar’s combined spectrum, 
as captured by Chandra and NuSTAR, can be represented by a 
simple power law that includes photoelectric absorption (see 
Figure 6). The resulting fitted parameters have been compiled in 
Table 3, with an Ny value of 2.31918 x 102? atoms cm~?. The 
fitting results with thawed Ny are consistent well with the results 
of Madsen et al. (2020), whose Ny was fixed to be 


spectra well, except x2, q of phase 0.3-0.4, which is 1.61, 
slightly big but still acceptable. There is no particular structure 
in the residual plot of this phase in Figure 7. 
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Figure 5. Panel (a): jointed spectrum around PSR J1811-1925 observed by Chandra and NuSTAR spectra in region (a) selected by NuSTAR. The gray-line indicates 
our spectral model: tbabs (powerlaw+vpshock); Panel (b): ratio between data and spectral model. 
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Figure 6. Panel (a): pulsed spectrum of PSR J1811—1925, combined Chandra and NuSTAR spectra. The dashed-line indicates a simple photoelectric absorbed power- 
law model. Panel (b): ratio between data and spectral model. 


We plot all photon index of phase-resolved spectra together From our results of phase-resolved spectra (Table 4 and Figure 
with pulse profile in Figure 8. The work of Ge et al. (2012) 8), the photon index I of PSR J1811-1925 shows no 
found that photon index of three young pulsars (PSR B0531 prominent evolution with phase, it keeps nearly constant at 
+21, PSR B1509-58, PSR B0540-69) evolves with phase. about 1.46 + 0.08, except the photon index of the off pulse 
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Table 3 
Fitting Results of Pulsar Spectrum 

No. Paame ë Vae 
1 Ny 2.317918 
2 Tt 1.33 + 0.06 
3 Norm (107^) 1.397913 
4 Flux” (1-10 keV) 0.91 + 0.04 
5 CFPMAB 1.672013 
6 Xa (d.o.f) 1.14(70) 
Notes. 


è Unit: 10°? atoms cm™°. 
> The absorbed X-ray flux unit is 107"? erg cm~? 1. 


Table 4 

Results of Phase-resolved Spectroscopy from NuSTAR Observation 
Phase li Creme Flux* (1-10 keV) ——-y2,4(d.o-f) 
0.0-0.1 1.45 +0.11 1.024813 0.38 0.03 0.92(74) 
0.1-0.2 1.38 0.09 1.017949 0.64 0.05 0.98(86) 
0.2-0.3 1.38 + 0.11 1.04+9-13 0.39+9:04 1.01(76) 
0.3-0.4 1.40 + 0.18 1.0370 72 0.22 + 0.03 1.28(57) 
0.4-0.9 1.81 + 0.08 1.007008 1.317908 1.03(186) 
0.9-1.0 1.547932 0.927938 0.16994 0.85(50) 
Notes. 


è The absorbed X-ray flux unit is 107}? erg cm~? s™!. 


phase, which is 1.81 + 0.08, softer than photon index of the 
pulse-on phase. The reason is that the pulse-off spectrum 
(phase 0.4-0.9) is dominated by the diffuse radiation of the 
PWN. The photon index of the pulse-off spectrum is also 
consistent with the PWN result of Roberts et al. (2003), whose 
PWN’s photon index is 1.54-1.90. 


4. Discussion and Summary 


The combined analysis of NuSTAR and Chandra data shows 
that the hard X-ray emission originates from PSR J1811-1925 
and its PWN. Its flux (1-10 keV) is 0.91 x 10°! erg cm~? s™!. 
Taking a distance of 5 kpc to the pulsar (Dean et al. 2008), its 
X-ray luminosity (1-10 keV) is 2.75 x 10% ergs’. This X-ray 
luminosity is similar to that of PSR B1509-58 (Gaensler et al. 
2002). Such pulsars with similar pulse profiles and energy 
spectra may belong to a new category, soft y-ray pulsars, as 
suggested by Kuiper & Hermsen (2015). Compared with y-ray 
pulsars, this type of pulsar has different properties: they all 
have a relatively young characteristic age (within tens of 
thousands of years); all have spin-down powers L,q above 
3 x 10°% ergs! (Kuiper & Hermsen 2015). 

However, within the soft y-ray pulsar population, they also 
have many different properties: 
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1. Pulse profile: PSR J1811-1925, PSR B1509-58, PSR 
J1846-0258, PSR J1930+1852, and PSR J1617—5055 
have a broad single peak. But Crab pulsar PSR B0531 
+21 has a double-peaked pulse profile, and Vela pulsar 
PSR B0833-45 even has a multiple-peaked pulse profile. 

2. Pulse phase dependence of the photon index: the spectra 
of PSR B1509-58 and PSR B0540-69 are harder at the 
peak of the pulse profile and softer at the wings. Although 
the spectra of PSR J1811—1925 have no apparent phase 
dependence considering the measured errors, we cannot 
rule out the possibility that its photon index evolves with 
the phase. Crab is different from these pulsars, which is 
softer at the two peaks (Ge et al. 2012). Cheng et al. 
(2000) proposed the CRZ model for outer magnetosphere 
of pulsars. In this model, pairs of discharges in the pulsar’s 
magnetosphere generate gamma-rays, followed by sec- 
ondary pairs that produce synchrotron radiation at the 
X-ray band. However, the view angle of PSR J1811-1925 
only allows us to detect the synchrotron radiation of the 
pairs created by the incoming particles. The outer gap 
model can explain the phase variation of the photon index 
of PSR B1509-58 (Hu et al. 2017). As shown in the results 
of Ge et al. (2012), where the photon index changes very 
little (within ~0.1) in different phases. While the spectra 
of PSR J1811-1925 show no apparent phase dependence, 
it is important to note that they have low statistics and 
large uncertainties. Therefore, we hope that NuSTAR can 
observe PSR J1811—1925 for an extended period of time. 

3. GeV radiation: PSR J1811-—1925, PSR J1846—0258, PSR 
J1930+1852, and PSR J1617-5055 have no GeV 
radiation detected (Wang et al. 2014), PSR B1509-58 
is not bright in GeV bands. Crab, Vela, PSR B0540-69 
and so on have significant GeV radiation. 


This class of sources (PSR J1811—1925, PSR B1509-58, PSR 
J1846-0258, PSR J1930+1852, and PSR J1617-5055) with 
similar X-ray energy spectra and pulse profiles, their photon 
index is around 1.2, their pulse profiles exhibit broad, 
structured single pulse. Therefore, such sources may be a 
new class of sources. 

Hard (21 keV) pulsed X-ray emission is believed to come 
from the magnetosphere of the pulsars powered by pulsar spin- 
down energy. A lot of models can explain this non-thermal 
pulsed X-/+-ray emission, and the most popular ones are outer 
gap model (Zhang & Cheng 1997; Hirotani 2015; Takata et al. 
2016) and polar gap model (Ruderman & Sutherland 1975). 
According to the outer gap model, it is suggested that the 
X-rays and soft y-rays emitted by pulsars may result from 
synchrotron radiation generated by electron pairs in the 
magnetic field. This type of radiation typically produces a 
wider pulse profile (Wang et al. 2013). In contrast, in the polar 
cap model, the X-rays and soft y-rays are thought to stem from 
curvature radiation, and the pulse profile tends to be narrower 
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Figure 8. Compare the relation between the photon index obtained in NuSTAR FPMA/FPMB and the pulse profile. The yellow region is the region of pulse off. 


(Wang et al. 2014). By our phase-resolved spectroscopy, we 
found the pulse radiation in a broad phase, which conflicts with 
the narrow phases of the curvature radiation. In addition, no 
GeV curvature photons detected suggest that curvature 
radiation may be outside our viewing angle (Wang et al. 
2014). Therefore, our results support that the pulsed radiation 
comes from synchrotron radiation. 

In summary, we have analyzed the data from the NuSTAR 
and Chandra observations of X-ray pulsar PSR J1811-1925. 
The following are our brief conclusions: 


1. From the X-ray pulse profile in different phases and 
regions, we can detect pulse radiation until 35-79 keV. 
The pulse profile does not show evolution as a function of 
energy. 
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2. Through its phase-resolved spectra, the photon index 
does not evolve as function of the pulse phase. For PSR 
B1509-58, no evidence of evolution was found at the 
beginning (Rots et al. 1998), but with more observations 
and the improvement of the signal-to-noise ratio, the 
photon index shows clear evolution (Ge et al. 2012). 
Therefore, we expect more observations of this source to 
obtain more accurate phase resolved spectra. 
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